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Abstract 



A search is presented for direct top squark pair production in final states with one isolated electron 
or muon, jets, and missing transverse momentum in proton-proton collisions at y/s = 7 TeV. The 
measurement is based on 4.7 fb _1 of data collected with the ATLAS detector at the LHC. Each top 
squark is assumed to decay to a top quark and the lightest supersymmetric particle (LSP). The data 
are found to be consistent with Standard Model expectations. Top squark masses between 230 GeV 
and 440 GeV are excluded with 95% confidence for massless LSPs, and top squark masses around 
400 GeV are excluded for LSP masses up to 125 GeV. 
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Weak scale Supersymmetry (SUSY) [l|-|9( is an exten- 
sion to the Standard Model (SM) that provides a solution 
to the hierarchy problem by introducing supersymmet- 
ric partners of all SM particles. In the framework of a 
generic i?-parity conserving minimal supersymmetric ex- 
tension of the SM (MSSM) [IM1, SUSY particles are 
produced in pairs, and the lightest supersymmetric par- 
ticle (LSP) is stable and can be a dark matter candidate. 
In a large variety of models, the LSP is the lightest neu- 
tralino, Xii which only interacts weakly and thus escapes 
detection. 

Light top squarks (stop) are suggested by naturalness 
arguments [3, E^- Searches for direct stop pair pro- 
duction have been previously reported by the CDF and 
DO experiments [121 EH- Searches for stops via ah pro- 
duction have been reported by the ATLAS [19h21] and 



CMS |22j, |23[ collaborations. In this Letter, one stop 



mass eigenstate (tj) is assumed to be significantly lighter 
than the other squarks. A search is presented for directly 
pair-produced stops, which are each assumed to decay to 
a top quark and the LSP. The final state for such a sig- 
nal is characterized by a top quark pair (ti) produced 
in association with large missing transverse momentum, 
the magnitude of which is referred to as £^ 1ISS , from the 
undetected LSPs. One top quark is assumed to decay 
hadronically and the other semileptonically. 

The ATLAS detector [24[ has a solenoid, surrounding 
the inner tracking detector (ID), a calorimeter, as well 
as a barrel and two endcap toroidal magnets support- 
ing the muon spectrometer. The ID consists of silicon 
pixel, silicon microstrip, and transition radiation detec- 
tors and provides precision tracking of charged particles 
for pseudorapidity \rj\ < 2.5 j25|. The calorimeter, placed 
outside the solenoid, covers \t}\ < 4.9 and is composed of 
sampling electromagnetic and hadronic calorimeters with 
either liquid argon (LAr) or scintillating tiles as the active 
media. The muon spectrometer surrounds the calorime- 
ters and consists of a system of precision tracking cham- 
bers in \t)\ < 2.7, and detectors for triggering in \rj\ < 2.4. 

The analysis is based on data recorded by the ATLAS 



detector in 2011 corresponding to 4.7 fb -1 of integrated 
luminosity with the LHC operating at a pp center-of-mass 
energy of 7 TeV. The data were collected using a logical 
OR of a E™ lss and a single-lepton (electron or muon) 
trigger. The combined triggers reach plateau efficiencies 
>98% for the chosen selection criteria on E™ lss and the 
leptons. Requirements that ensure the quality of beam 
conditions, detector performance and data are imposed. 

Monte Carlo (MC) event samples using the full AT- 
LAS detector simulation [26| based on the Geant4 pro- 
gram [2jj are used to aid in the description of the 
background and to model the SUSY signal. The ef- 
fect of multiple pp interactions per bunch crossing is 
also simulated [28}. Production of top quark pairs is 
simulated with MCSNL0 4.01 [H, [H, alternatively us- 
ing ALPGEN 2.14 [31] and PowHeg HVQ patch 4 [32- 
|34| . The data modeling is improved for high jet mul- 
tiplicities by reweighting the MCQNL0 sample to match 
the jet multiplicity distribution in ALPGEN. Uncertain- 
ties associated with initial- and final-state radiation 
(ISR/FSR) [H] are assessed using AcerMC 3 . 7 [II] sam- 
ples. A top quark mass of 172.5 GeV is used consis- 
tently. W and Z/"f* production in association with 
jets are each modeled with ALPGEN. Diboson VV (WW, 
WZ and ZZ) production is simulated with ALPGEN and 
cross-checked with HERWIG 6.520 (37|. Single top pro- 
duction is modeled with MC0NL0, and ti events pro- 
duced in association with Z, W or WW (ti + V) 



are generated with MADGRAPH 5 38]. Next-to- leading- 
order (NLO) PDFs CT10 [!§] are used with all NLO 
MC samples. For all other samples, LO PDFs are 
used: MRSTmcal with HERWIG, and CTEQ6L1 [Il| with 
ALPGEN and MADGRAPH. Fragmentation and hadronization 
for the ALPGEN and MC0NL0 samples are performed with 
HERWIG, using JIMMY 4.31 [12] for the underlying event, 
and for the MADGRAPH samples PYTHIA 6 . 425 [43[ is used. 
The ti, single top and ti + V production cross sections 
are normalized to approximate next-to-next-to-leading 
order (NNLO) [4J], next-to- next-to-leading - log arithmic 
accuracy (NLO+NNLL) and NLO 0] calcula- 
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tions, respectively. QCD NNLO FEWZ [Ig| inclusive W 
and Z cross sections are used for the normalization of 
the W+jets and Z+jets processes. Expected diboson 
yields are normalized using NLO QCD predictions ob- 
tained with MCFM (50ll5H. 

Stop pair production is modeled using 
Herwig++ 2.5.2 [H^]. The l Y is mostly the part- 
ner of the right-handed top quark, and the Xi is almost 
a pure bino. A signal grid is generated with a step 
size of 50 GeV both for the stop and LSP mass values. 
Signal cross sections are calculated to NLO in the strong 
coupling constant, including the resummation of soft 
gluon emission at next-to-leading-logarithmic accuracy 
(NLO+NLL) [5314551 ] . The nominal cross section and the 
uncertainty are taken from an envelope of cross section 
predictions using different PDF sets and factorization 
and renormalization scales [56]. The t 1 t 1 cross section 
for m ?i = 400 GeV is (0.21 ± 0.03) pb. 

Events must pass basic quality criteria to reject detec- 
tor noise and non-collision backgrounds [57, 58 1 and are 
required to have > 1 reconstructed primary vertex asso- 
ciated with five or more tracks with transverse momen- 
tum pt > 0.4 GeV. Events are retained if they contain 
exactly one muon [H^| with |??|_< 2.4 and pt > 20 GeV 
or one electron passing 'tight' [60j selection criteria with 
\r]\ < 2.47 and px > 25 GeV. Leptons are required to 
be isolated from other particles. The scalar sum of the 
transverse momenta of tracks above 1 GeV within a cone 
of size AR = 0.2 around the lepton candidate is required 
to be < 10% of the electron pr, and < 1.8 GeV for the 
muon. Events are rejected if they contain additional lep- 
tons passing looser selection criteria [6ll ] . Jets are recon- 
structed from three-dimensional calorimeter e nerg y clus- 
ters using the anti-fc t jet clustering algorithm 62] with a 
radius parameter of 0.4. The jet energy is corrected for 
the effects of calorimeter non-compensation and inhomo- 
geneities using pt- and ry-dependent calibration factors 
based on MC simulations and validated with extensive 
test-beam and collision-data studies [63| ■ To suppress jet 
background originating from uncorrelated soft collisions, 
> 75% of the summed pt of all tracks associated to a 
jet must come from tracks associated to the selected pri- 
mary vertex. Events with four or more jets with \rj\ < 2.5 
and pt > 80, 60, 40, and 25 GeV are selected. At least 
one jet needs to be identified as a &-jet, which is a jet 
containing a 6-hadron decay. These are identified using 
the 'MV1' 6-tagging algorithm [64J which exploits both 
impact parameter and secondary vertex information. An 
operating point is employed corresponding to an average 
75% 6-tagging efficiency and to a < 2% misidentification 
rate for light-quark/gluon jets for jets with pt > 20 GeV 
and |?7| < 2.5 in ti MC events. 

Ambiguities between overlapping leptons and jets are 
resolved by discarding either the jet or lepton candi- 
dates 0] based on their separation AR. The measure- 
ment of E!f lss is based on the transverse momenta of all 



electron and muon candidates, all jets after overlap re- 
moval, and all calorimeter energy clusters not associated 
to such objects. The background is reduced by requiring 
A(j) m i n > 0.8, where A</> m ; n is the minimum azimuthal 
separation between the two highest pt jets and the miss- 
ing transverse momentum direction. A requirement on 
the three-jet mass rrijjj of the hadronically decaying top 
quark specifically rejects the dileptonic ti background, 
where both W bosons from the top quarks decay lepton- 
ically. The jet-jet pair having invariant mass > 60 GeV 
and the smallest AR is selected to form the hadronically 
decaying W boson. The mass rrijjj is reconstructed in- 
cluding a third jet closest in AR to the hadronic W boson 
momentum vector and 130 GeV < rrijjj < 205 GeV is re- 
quired. 

Five signal regions (SR A-E) are defined in order to op- 
timize the sensitivity for different stop and LSP masses. 
For increasing stop mass and increasing mass difference 
between stop and LSP the requirements are tightened on 
Ef iss , on the ratio Ef iss / \fWr~ ', where H T is the scalar 
sum of the momenta of the four selected jets with high- 
est pt, and on the transverse mass rax 65J, as shown in 
Table fl] The number of observed events in each SR after 
applying all selection criteria are given in Table HIl 

TABLE I. Selection requirements defining the SR A-E. 
Requirement SR A SR B SR C SR D SR E ~ 



r l/2- 



m T [GeV] > 



150 


150 


150 


225 


275 


> 7 


9 


11 


11 


11 


120 


120 


120 


130 


140 



The product of the kinematic acceptance, detector and 
reconstruction efficiency (A ■ e) varies between 4% and 
0.3% for SR A and between 3% and 0.01% for SR E as 
the stop-LSP mass difference varies between 550 GeV and 
250 GeV. 

The dominant background arises from dileptonic ti 
events in which one of the leptons is either not identi- 
fied, is outside the detector acceptance, or is a hadron- 
ically decaying r lepton. In all these cases, the ti de- 
cay products include two or more high-px neutrinos, re- 
sulting in large E™ ss and m-p. Three control regions 
(CRs) enriched in dileptonic ti events (2-lep TR), single- 
leptonic ti events (1-lep TR), and VF+jets events (1-lep 
WR) are designed to normalize the corresponding back- 
grounds using data. The 2-lep TR differs from the SRs 
by selecting events with exactly two leptons, applying no 
requirements on wit, E T 11SS / \J Ht and rrijjj, and by re- 
quiring E^ iss > 125 GeV. The 1-lep TR and 1-lep WR 
have selection criteria identical to SR A, except the mx 
requirement is changed to 60 < rax < 90 GeV and the 1- 
lep WR has a 6-jet veto instead of a 6-jet requirement, ti 
production accounts for > 90% of events in the top CRs 
and VL+jets production for > 60% in the W CR. The 
signal contamination reaches a maximum of 8% in the 
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2-lep TR for mj = 200 GeV. The multijet background, 
which mainly originates from jets misidentified as lep- 
tons, is estimated using the matrix method [6l| . Other 
background contributions (VV, tt + V and single top) are 
estimated using MC simulation normalized to the theo- 
retical cross sections. The Z+jets background is found 
to be negligible. 

Good agreement is observed between data and the SM 
prediction before using the CRs to normalize the ti and 
VF+jets backgrounds. As an example, Fig. Q] shows the 
agreement of the £"™ lss distributions in the 2-lep TR, 
and the E™ lss distribution in SR A. In addition, the hit 
distribution for a looser requirements region — E^ lss > 
40 GeV and no requirements on E™ lss / \/Ht and rrijjj 
(pre-selection) — is shown. 

Simultaneous fits to the numbers of observed events in 
the three CRs and one SR at a time are performed to nor- 
malize the ti and W+jets background estimates as well 
as to search for an excess from a potential signal contri- 
bution. The 1-lep and 2-lep TRs have ti normalizations 
that float independently and that are found to be in good 
agreement with each other. The ti estimates in the SRs 
are based on the 2-lep TR, as this minimizes the extrap- 
olation uncertainties in the fit. Systematic uncertainties 
are treated as nuisance parameters with Gaussian prob- 
ability density functions. 

The dominant systematic uncertainties in the fitted ti 
background estimate are theoretical and modeling uncer- 
tainties. They are determined by using different gener- 
ators (MCONLO, PowHeg and ALPGEN), different showering 
models (HERWIG and PYTHIA) and by varying ISR/FSR 
parameters, and amount to 10-30% on the extrapolation 
from the CR to the various SRs. Electroweak single top 
production is associated with an 8% theoretical uncer- 
tainty I45l-l47l ] and the ti+V background has a 30% uncer- 
tainty [ig. The difference between ALPGEN and HERWIG 
predictions is used to assess the uncertainty on the di- 
boson background, and the uncertainty on the multijet 
background is based on the matrix method. Both of these 
uncertainties are estimated as 100%. 

Experimental uncertainties affect the signal and back- 
ground yields estimated from MC events and are dom- 
inated by uncertainties in the jet energy scale, jet en- 
ergy resolution and 6-tagging. Uncertainties related to 
the trigger and lepton reconstruction and identification 
(momentum and energy scales, resolutions and efficien- 
cies) give smaller contributions. Other small uncertain- 
ties are due to modeling of multiple pp interactions, the 
integrated luminosity, and the limited numbers of MC 
and data events. The uncertainty on A ■ e varies between 
9% and 16% as the simulated stop-LSP mass difference 
varies between 550 GeV (SR E) and 250 GeV (SR A and 
B). 

Table UU shows the results of the background fit to the 
CRs, extrapolated to the SRs. The fitted numbers of ti 
and W+jets events are compatible with the MC predic- 
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FIG. 1. Top: 
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distributions for the 2-lep TR. Center: 



the pre-selection requirements (see text). All plots show the 
electron and muon combined channel, before normalization 
fits. Hatched areas indicate the combined uncertainty due to 
MC sample size and the jet energy scale. 



tions, with factors of 1.01 and 0.90 applied, respectively. 
To assess the agreement between the SM expectation and 
the observation in the SRs, a second set of simultaneous 
fits including one SR at a time and all CRs is performed. 
The pa-valnes (probing the background-only hypothesis) 
obtained are given in Table [TTJ No significant excess of 
events is found. 

One-sided exclusion limits are derived using the CL S 
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TABLE II. Numbers of observed events in the five signal regions and three background control regions, as well as their estimated 
values and all (statistical and systematic) uncertainties from a fit to the control regions only, for the electron and muon combined 
channel. The expected numbers of signal events for = 400 GeV (500 GeV) and m^o = 1 GeV for benchmark points 1 (2) 
are listed for comparison. The central values of the fitted sum of backgrounds in the control regions agree with the observations 
by construction. Furthermore, po-values and 95% CL S observed (expected) upper limits on beyond-SM events are given, using 
simultaneous fits including one SR at a time and all CRs. 



Regions 


SR A 


SR B 


SR C 


SR D 


SR E 


2-lep TR 


1-lep TR 


1-lep WR 


it 

tt + V , single top 
V+jets, VV 
Multijet 


36 ±5 
2.9 ±0.7 
2.5 ± 1.3 
0.4 ±0.4 


27 ±4 
2.5 ±0.6 
1.7 ±0.8 
0.3 ±0.3 


11 ± 2 
1.6 ±0.3 
0.4 ±0.1 
0.3 ±0.3 


4.9 ± 1.3 
0.9 ±0.3 
0.3 ±0.1 
0.3 ±0.3 


1.3 ±0.6 
0.4 ±0.1 
0.1 ±0.1 

o o+ - 3 

u - u -o.o 


109 ± 10 

7.2 ± 1.3 

1.6 ± 0.8 

0+ ' 6 
u - u -o.o 


364 ± 23 
18 ±3 
38 ±11 

1.7 ±1.7 


59 ± 19 
6.1 ± 1.6 
162 ± 23 
0.8 ±0.8 


Total background 
Signal benchmark 1 (2) 
Observed events 


42 ±6 
25.6 (8.8) 
38 


31 ±4 
23.0 (8.1) 
25 


13 ± 2 
17.5 (6.9) 
15 


6.4 ± 1.4 
13.5 (6.2) 
8 


1.8 ±0.7 
7.1 (4.5) 
5 


118 ± 10 
1.7 (0.6) 
118 


421 ± 20 
2.3 (0.6) 
421 


228 ± 15 
0.4 (0.1) 
228 



po-values 0.50 0.50 0.32 0.24 0.015 

Obs. (cxp.) JVbsm < 15.1 (17.2) 10.1 (13.8) 10.8 (9.2) 8.4 (7.0) 8.2 (4.6) 



method [66| , based on the same simultaneous fit method 
but taking the predicted signal contamination in the CRs 
into account. To obtain the best expected combined ex- 
clusion limit, a mapping in the stop-LSP mass plane is 
constructed by selecting the SR with the lowest expected 
CL S value for each grid point. The expected and ob- 
served 95% CL S exclusion limits are displayed in Fig. [2j 
Stop masses are excluded between 230 GeV and 440 GeV 
for massless LSPs, and stop masses around 400 GeV are 
excluded for LSP masses up to 125 GeV. These values 
are derived from the — 1 of observed limit contour. 
These stop mass limits significantly extend previous re- 
sults. Limits on beyond-SM contributions are derived 
from the same simultaneous fit but without signal model- 
dependent inputs (i.e. without signal contamination in 
the CRs, and without signal systematic uncertainties). 
The resulting limits are shown at the bottom of Table ILT1 
In summary, a search for stop pair production is pre- 
sented in final states with one isolated lepton, jets, and 
missing transverse momentum in ^fs = 7 TeV pp colli- 
sions corresponding to 4.7 fb -1 of ATLAS 2011 data. 
Each stop is assumed to decay to a top quark and a 
long-lived undetected neutral particle. No significant ex- 
cess of events above the rate predicted by the Standard 
Model is observed and 95% CL S upper limits are set on 
the stop mass in the stop-LSP mass plane, significantly 
extending previous stop-mass limits. 



ACKNOWLEDGEMENTS 



^ 300 
CD 

^250 



200 



Expected limit (±1 o e . 



ATLAS 

1-lepton + jets + E' 

J L dt = 4.7 fb Js=7 TeV .'0.72 0.16 0.12 0.10 



150 



100 



50 



production, t t -» t x 

n 1 i 1 



"~ 1 i 1 i 1 r 

Observed limit (±1 of 



All limits at 95% CU 




200 250 300 350 400 450 500 550 600 

ny [GeV] 



FIG. 2. Expected (dashed) and observed (solid curve) 95% 
CL S excluded region (under the curve) in the plane of m-o 
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limit on the signal cross section, in pb. 
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